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A renormalization group method is developed with which 
thermodynamic properties of a weakly interacting, confined 
Bose gas can be investigated. Thereby effects originating 
from a confining potential are taken into account by periodic 
boundary conditions and by treating the resulting discrete en- 
ergy levels of the confined degrees of freedom properly. The 
resulting density of states modifies the flow equations of the 
renormalization group in momentum space. It is shown that 
as soon as the characteristic length of confinement becomes 
comparable to the thermal wave length of a weakly inter- 
acting and trapped Bose gas its thermodynamic properties 
are changed significantly. This is exemplified by investigat- 
ing characteristic bunching properties of the interacting Bose 
gas which manifest themselves in the second order coherence 
factor. 



I. INTRODUCTION 

The experimental realization of Bose Einstein conden- 
sation of trapped, ultra cold, weakly interacting atomic 
gases has renewed the interest in their thermody- 
namic properties. Definitely, many of these properties 
have already been studied more than thirty years ago 
in connection with the theoretical efforts to describe su- 
perfluidity of strongly interacting helium atoms. How- 
ever, in these recent experimental realizations the physics 
of the trapping process and of the internal structure of 
the trapped bosons still pose interesting new theoretical 
questions which are not yet understood completely 

So far the majority of theoretical approaches aiming at 
describing thermodynamic properties of weakly interact- 
ing, trapped Bose gases has concentrated on mean field 
approaches jjj. At zero temperature these mean field 
approaches lead to the well known Gross-Pitaevski equa- 
tion and to their multicomponent generalizations. These 
theoretical descriptions have proven useful in obtaining 
first quantitative understandings of many recent trapping 
experiments [El. However, it is known that mean field 
theories do not yield an accurate description of thermo- 
dynamic properties of interacting gases close to the tran- 
sition point of a phase transition Furthermore, with 
improving experimental accuracy it is expected that in 
the near future possible deviations from mean field results 
might become accessible experimentally even in regimes 
far away from transition points. 



Motivated by these prospects in this paper a theo- 
retical description of thermodynamic properties of con- 
fined, weakly interacting Bose gases is developed which 
is based on renormalization group (RG) methods Hf). 
RG methods are known to yield accurate descriptions of 
thermodynamic partition functions also in regions close 
to a second order phase transition. Furthermore, by 
these methods certain universal quantities, such as crit- 
ical exponents, might even be evaluated by analytical 
means. Typically the evaluation of thermodynamic par- 
tition functions with the help of RG methods is only 
slightly more complicated than by mean field methods. 
Thus, these RG methods represent a convenient theoreti- 
cal approach complementing elaborate numerical Monte- 
Carlo simulations of partition functions with which one 
may also test limitations of mean field theories. 

In the special case of an unconfined, homogeneous 
three-dimensional interacting Bose gas such a RG ap- 
proach has already been developed recently by Bijlsma 
and Stoof J8|. However, their quantitative predictions 
of the dependence of the transition temperature on the 
scattering length, for example, have been more than one 
order of magnitude larger than corresponding numerical 
predictions which are based on elaborate Monte-Carlo 
simulations of the relevant partition function |^]. In a 
recent perturbative approach this discrepancy has been 
clarified |Hj|nj]. These perturbative studies are in sat- 
isfactory agreement with the RG calculations of Bijlsma 
and Stoof J8| and with recent experiments |l2| and at- 
tribute the numerical inaccuracies of the previous Monte- 
Carlo approach to incorrect extrapolation procedures. 
This example demonstrates that despite extensive pre- 
vious work there are still open questions even in the con- 
text of the traditional case of a homogeneous and un- 
confined interacting Bose gas |p^| . In addition, the ex- 
perimental possibility to control and manipulate trapped 
Bose condensates has raised many new interesting phys- 
ical questions most of which have not been tackled yet 
theoretically beyond the framework of mean field approx- 
imations. 

In this paper a RG approach is developed which 
is capable of determining thermodynamic properties of 
weakly interacting and confined Bose gases. This ap- 
proach aims at describing physical effects of spatial con- 
finement beyond the local density approximation by 
modeling effects originating from confinement by periodic 
boundary conditions. However, contrary to the usual 
continuum-limit-approximation within our approach the 
required summations over all possible states of the atomic 
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center of mass motion are not replaced by integrations 
but the discrete energy levels resulting from the confined 
degrees of freedom are taken into account. Thus the pre- 
sented approach should yield a valid description even in 
cases in which the thermal wave length of an interacting 
Bose gas becomes comparable to the characteristic length 
of confinement and in which the local density approxima- 
tion breaks down. Indeed in the subsequent treatment it 
is demonstrated that as soon as the confinement length 
becomes comparable to the thermal wave length of the in- 
teracting Bose gas thermodynamic properties of an inter- 
acting Bose gas are changed significantly. This character- 
istic behavior is exemplified by investigating various ther- 
modynamic properties of an interacting Bose gas, such as 
its bunching properties which manifest themselves in the 
second order coherence factor and which are accessible to 
experimental observation. Though the presented RG ap- 
proach differs slightly from the one developed by Bijlsma 
and Stoof § previously it is demonstrated that resulting 
critical properties, such as the dependence of the tran- 
sition temperature on the scattering length, are in sat- 
isfactory agreement with these previous RG results and 
with the recent perturbative calculations of Holzmann et 
al. [[00. 

The paper is organized as follows: In Sec. II the gen- 
eral theoretical approach is developed. Starting from 
the path integral representation of the grand canonical 
partition function an RG method is developed for its 
evaluation. As effects of a confining potential are taken 
into account by periodic boundary conditions which still 
guarantee translational invariance of the problem the RG 
procedure is developed in momentum space. As a main 
approximation the random phase approximation (RPA) 
|L4j is used which is applied to the symmetry broken 
phase. In Sec. Ill resulting numerical results are pre- 
sented. Firstly, we concentrate on characteristic prop- 
erties of an unconfined homogeneous three-dimensional 
Bose gas which have been of interest recently. Secondly, 
the influence of confinement is exemplified by investigat- 
ing isothermal properties of the pressure and of the sec- 
ond order coherence factor of an interacting Bose gas. 

II. THEORETICAL TREATMENT 

In this section a RG approach is developed for the eval- 
uation of the grand canonical partition function of an in- 
teracting and confined Bose gas in D spatial dimensions. 
Effects of confinement are taken into account by periodic 
boundary conditions and by proper summations over the 
discrete energies associated with the confined degrees of 
freedom. Accordingly, the RG approach is developed in 
momentum space. The RG flow equations are derived for 
the temperature regime below the phase transition with 
the help of the random phase approximation (RPA) . 
This dynamical regime is characterized by a non-zero or- 
der parameter. Details of the derivation of these RG 



equations are given in appendix A. For the sake of com- 
pleteness the RG equations which apply to the dynamical 
region above the phase transition, i.e. in the region of a 
vanishing order parameter, are summarized in appendix 
B. These latter equations been already discussed in detail 
previously |qjl4|]. 

A. Path integral representation of the partition 
function 

Thermodynamic properties of an interacting system of 
bosons can be described by its (grand-canonical) parti- 
tion function (T^] 

Z(z, 0) = Tt[ e -m-^ (i) 

with the fugacity z = (n denotes the chemical po- 
tential) and with the inverse temperature f3 = l/{kT) 
(k and T are the Boltzmann constant and temperature). 
The Hamiltonian and the number operator of the inter- 
acting bosons are denoted H and N, respectively. For 
interacting bosons this partition function can be repre- 
sented as a functional integral over a complex scalar field 
<I>. In D spatial dimensions this path integral represen- 
tation is given by |l5|,f7j 

Z(z, (3) = J dm, 0), er]e- s (*<**> (2) 

with the (dimensionless) action 

£($,$*)=/ / ded D ^{$*(e,e)h-A 5 - 

JO JVA D Z 

\ ( d D e\m,e)\ 2 [PAV(^-^)}\<f^',e)\ 2 }. (3) 

^ JVA D A 

The volume within which the Bose system is confined 
is denoted V. For the sake of simplicity this volume 
will be assumed to be of cubic shape and the resulting 
spatial confinement will be described by imposing peri- 
odic boundary conditions on the complex field $(£, 0). 
In Eq.(||) scaled quantities have been introduced which 
involve a yet arbitrary characteristic momentum (HA) 
whose associated inverse temperature is given by (3\ = 
m/(h 2 A 2 ). The mass of the interacting bosons is de- 
noted to. The scaled chemical potential M is defined 
by M — fi(3\. The two-body interaction potential be- 
tween the bosons is denoted ^(x — x') and £ = xA 
and = T/(Tif3\) are scaled spatial coordinates and 
the scaled imaginary time. The fugacity z appearing in 
Eq.(||) is related to the chemical potential M by 

In* = M - ~/?aV(x = 0). (4) 

The integration measure appearing in Eq.(|^) is defined 
later (cf. Eq.(|)). 
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Due to the trace-operation involved in the evaluation 
of the partition function of Eq.(Q) the functional inte- 
gration in Eq.(||) has to be performed over all complex 
fields <&(£, 9) which are periodic in the imaginary scaled 
time 9, i.e. $(£,6) = $(£,9 + 0//3 A ) §■ If one also 
imposes spatial periodic boundary conditions which still 
preserve the translational invariance of the problem it is 
convenient to transform the functional integral of Eq. (j|) 
into the momentum - frequency representation. This is 
achieved by Fourier transforming the complex fields ac- 
cording to 



*(f,e) = X;[ 



][- 



D — iuj n 



\9V: 



(5) 



with the Matsubara frequencies u„ = 2irnf3 A / f3 and with 

the scaled wave vectors km = 27rm.W / (i»A) ( mW being 
integer). Thereby the quantities Li (i = 1, ...,D) denote 
the lengths of the confining cubic volume. Defining (hA) 
as the ultraviolet momentum cut-off of the path integral, 
which is much larger than all physically significant mo- 
menta, implies that the scaled wave vectors appearing in 
Eq.(||) are restricted to the region < |k m | < 1. Insert- 
ing this Fourier decomposition into Eq.(^) one obtains 

**) - 5}-tw„ + tm - M}\ip mn \ 2 + (6) 

mil 

(VAVy^p/pA)- 1 - G(ki+k 3 )(^_ mi _ ni r x 

mini ;. . . ;rri4n4 

( l P- m2 -n 2 )*<P m3 n 3 <Pn U nJ{mi + m 2 + m 3 + m 4 ) X 

5(m +ri2 + TI3 + ru) (7) 

with e m = (S denotes the Kronecker-delta func- 

tion.) The two-particle interaction is characterized by its 
Fourier components 



G(k) = / d D ^(3 A V(^/A). 

JVA D 



(8) 



In the subsequent discussion we assume that the two 
particle interaction is of short range so that G(k) is inde- 
pendent of momentum, i.e. G(k) = G. Thus, introducing 
the momentum cut-off (hA) the path integral of Eq 
becomes R.O 



Z(z,p) = [H J 



The normalization factors 



(9) 



(10) 



with 



2sinh(C m /2) = e 



_ JP/0 A )(e m -M)/2 



guarantee that in the limit of vanishing interactions, i.e. 
G = 0, the partition function reduces to the well known 
expression for the ideal Bose gas M, i.e. 



Z(z, (3) = 1[{[1 



-(/VAvXem-M)! 



(11) 



As expected on physical grounds, for fl/(3 A ,l/M ^> 1 
Z(z, (3) should become independent of the momentum 
cut-off (HA). 9bviously this requirement is fulfilled for 
the ideal gas. 



B. Renormalization group approach for the 
evaluation of the partition function 

One of the simplest methods of evaluating the parti- 
tion function of Eq.Q is the mean field approximation 
]l7|] . Thereby one expands the field 9) appearing in 
Eq.(||) quadratically around the most probable, uniform 
field configuration 



<p=y/(VA D )(0/l3A)M/G 



(12) 



which is determined by the requirement 

SS &**)\*=*.=?, y / (VA°)(p/ M = °- The ^suiting 
Gaussian integrals of this quadratic expansion can be 
performed in a straightforward way. For G > this most 
probable field configuration is nonzero only if M > 0. 

A more sophisticated and more accurate method of 
evaluating the partition function makes use of renormal- 
ization group methods. The basic idea of the renor- 
malization group procedure is to perform the integra- 
tions in Eq.(^|) successively ||. In each step only field 
components Stp mn are integrated out whose momenta 
are located in an infinitcsimally small momentum shell 
around the maximum momentum (HA), i.e. for which 
< |k m | < 1 with < I -C 1. All other small- 
momentum field components <p mn outside this momen- 
tum shell, which constitute the field <p < , are left un- 
changed. If only an infinitesimal momentum shell is 
integrated out, a quadratic expansion of S , ($, <£*) with 
respect to the fields Sip mn is sufficient for the further 
evaluation of the partition function jl6| . 

The desired aim of this integration over the large- 
momentum field components is to obtain a new scaled 
partition function which is similar to the original one ex- 
cept for possible renormalizations of the characteristic 
parameters M and G ■ I n particular this new scaled 
partition function has to have the same momentum cut- 
off (hA). By repeated application of this transformation 
one eventually integrates out all momentum components 
and obtains the value of the partition function. However, 
in general this desired aim can be achieved only approxi- 
mately. Two commonly used approximation schemes in- 
volved are perturbation theory |5[ and the random phase 
approximation (RPA) ]l4|] . In our subsequent treatment 
the RPA is employed. It is particularly useful in the case 
of a non- vanishing most probable configuration Tp |^| . 

In the subsequent development we separate the zero- 
temperature contribution 



3 



from the grand thermodynamic potential according to 

n(M,(3)=-\nZ(z,P)/(p/(3 A ) 



n(M, (3) = n(M, (3^oo)+ w(M, (3). 



(13) 



For the zero-temperature contribution f2(Af, (3 — > oo) ac- 
curate approximations are available. In the simplest form 
of the mean field approximation this zero-temperature 
contribution is given by 17 



fl(M,/3 



oo 



(14) 



Note that the depletion term |Q is not included in 
Eq. (O) . Typically this term is negligibly small for a 
weakly interacting Bose gas. 

It is shown in appendix A that on integrating out mo- 
mentum components of the fields tp mn which are located 
in the momentum shell e~ l < |k m | < 1 (0 < I <C 1) the 
RPA yields the differential scaling relation 



du(l) (VA U ) 



d{l)e~ m \n{l - e- x{1 ^) 



(15) 



oo). Thereby the quantity 



dl (f3/f3 A ) 

with w(M,0) = u(J, - 
A(Z) = (/3(/)//3 A )Ve>(e> +2M(/)) reflects the Bogoli- 
ubov dispersion relation at each step of the renormaliza- 
tion with the scaled inverse temperature (3(1) = (3e~ 21 . 
The differential scaling of the chemical potential M(l) 
and of the coupling strength G(l) is governed by the RG 
equations 



dM(l) 
dl 



2M(l) + 

Jcoth(A(Q/2) - 1] 
d{l)((3(l)/f3 A ) x 

\2M(l)G(l) - G(/)2e>] - 
3 M(l)G(l) ^ 



d{l){(3{l)/(3 A 
1 



{ 



2sinti (A(Z)/2) 
coth(A(Z)/2) - 1] 

W) 



2A(0 2 
+ 



}[2e >+ M(Z)] 5 



(16) 



dGQ) 
dl 



-(D-2)G(l) 
d(l)([3(l)/[3 A 



[coth(A(0/2) - 1] 



d(l)(P(l)/(3 A f 



G{lf 



2A(0 



3G(Z) 2 



2X(l) 2 L 2sinh 2 (A(0/2) 
coth(A(0/2) - H 



A(Z) 



L}[2c> + M(0^ 



(17) 



The (scaled) energy of the eliminated (infinitesimal) mo- 
mentum shell is denoted e> = 1/2. The number of states 



per unit volume inside such an infinitesimal momentum 
shell e~ l < |k m | < 1 (0 < I < 1) is given by [d{l){A D l)}. 
In the continuum limit, in D spatial dimensions one ob- 
tains the I— independent density 



d(l) 



n, 



(2ir 



if 



(18) 



M(l 


-0) 


= M, 


G(l 


= 0) 


= G, 


uj(1 


= 0) 


= 0. 



with n D = 2tt d / 2 /T(D/2) denoting the surface of a D- 
dimensional sphere, (r denotes the Euler-Gamma func- 
tion.) 

Eqs.(|l5|- |l^) constitute the set of differential RG equa- 
tions from which the temperature dependent part of the 
grand thermodynamic potential lj(M, (3) = uj(l — > oo) 
can be evaluated. They have to be solved subject to the 
initial conditions 



(19) 



The RG equations (15-17) are valid for an arbitrary 
number of spatial dimensions D of the interacting Bose 
gas as long as effects of confinement can be described 
by periodic boundary conditions. Effects of confinement 
are characterized by the quantity d(l) which reduces to 
the result of Eq.(|l§|) in the continuum limit, i.e. in the 
absence of any confining influence. If one or more degrees 
of freedom are confined, d(l) is modified as one has to take 
into account the resulting discrete energy levels in the 
confined degree of freedom. In a case where one degree 
of freedom of the interacting Bose gas is confined, for 
example, one obtains 



d(l) 



zMAe 



(^F2 2. J ^ 

y 1 M=-oo J ~ L 

fl D Tr + 2n[L z e- l A/(2n)} 



{2ir) D 



L 7 e- l K 



(20) 



Thereby L z is the length of the confining (square-well) 
potential and [x] denotes the largest integer which is 
smaller or equal to x. 

For dimensions less than or equal to two one has to 
reconsider the validity of the k-independent approxima- 
tion G(k) = G of the interparticle potential. If one 
assumes a spherically symmetric short-ranged interpar- 
ticle interaction then at low energies one may expand 
G(k) = G + G 2 k 2 + G 4 k 4 + .... It is straightforward 
to demonstrate that for / — * oo (which implies (3(1) = 
f3e~ 21 — * 0) the trivial scaling of Gk (k = 0,2,4,...) is 
given by G k (l) = G fc e[ 2 ~ (D ~ 2+fc)li • This implies that al- 
ready at D = 2 not only Go but also G2 is no longer 
irrelevant. Thus for a reliable description of the temper- 
ature dependence of two-dimensional Bose gases also the 
k-dependence of the interparticle interaction has to be 
taken into account at least to lowest order. 
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It is worth mentioning that Eqs.(|l5|- |l7]) reduce to the 
corresponding mean field results if one neglects fluctua- 
tions of the most probable configuration Tp. Formally this 
is equivalent to setting d(l) = (or to putting A(<£><) = 
in Eq.(|34|) of appendix A). In this case the scaling of the 
characteristic parameters M(l) and G(l) is governed by 
the trivial scaling which depends on the dimension D 
of the system only. Furthermore, these mean field re- 
sults also appear in the limit of zero temperature where 

All thermodynamic properties which are derivable 
from the grand thermodynamic potential f2(M, /?) can 
be determined by these RG equations. Thus pressure P 
and number of particles N, for example, are determined 
by the relations 



kT 
V 



MM,f3)^(2n(3/M D /\ (21) 
BO 

-(2ttP/(3 a ) d ^ — (M,(3) (22) 



with the thermal wave length \th 



2irh 2 /(mkT) 



y/2Trl3//3\/A. Another physical quantity which is of cur- 
rent experimental interest is the spatially averaged sec- 
ond order coherence factor of the interacting Bose gas. 
In terms of the conjugate quantized field operators f/>(x) 
and ifi' (x) of the Bose gas it is defined by 



<7 (2) (0) 



^Jd D x < ^t(x)^t(x)y3(x)yi(x) > 
[i fd D x< V>t(xMx) >] 2 



(23) 



with < . > denoting averaging over the thermodynamic 
equilibrium distribution of the interacting Bose gas. This 
definition is completely analogous to the corresponding 
quantity in the quantum optical context. As apparent 
from Eq.(l) this quantity is related to the grand thermo- 
dynamic potential of Eq.(14) by 



gW(0) = 2(VA D ) §^liS. , 
[ ' 1 j [|« ( M,/?,G)]2 



(24) 



In the partial derivatives entering Eq. (|2J) the grand 
thermodynamic potential f2 has to be considered as a 
function of the fundamental parameters (M, j3, G). In 
order to keep the notation as simple as possible this ex- 
plicit dependence on G has not been indicated explicitly 
in all other previous equations. 

Besides taking into account effects originating from 
confinement the RG flow equations (15-17) differ from 
the ones derived previously by Bijlsma and Stoof Q also 
in two other respects. First of all, separating off the 
zero-temperature contribution from the grand thermody- 
namic potential according to Eq.(13) changes the prop- 
erties of the flow equations. These changes are partic- 
ularly prominent in the region where X(l) 3> 1, i.e. in 
the regions where the value of the effective momentum 
cut-off is still large. The RG equations (15-17) have the 



property that for X(l) ^> 1 the influence of the nontrivial 
scaling which is proportional to the characteristic level 
density d(l) is exponentially small. Secondly, evaluating 
the particle density directly from the grand thermody- 
namic potential according to Eq. ( p2| ) yields results which 
are slightly different from the ones obtainable from the 
approximate flow equation used in Ref. S . 



III. NUMERICAL RESULTS 

In this section thermodynamic properties of a con- 
fined interacting Bose gas are investigated on the basis 
of the RG approach of Sec. II. We concentrate on cases 
in which one degree of freedom of a three-dimensional 
Bose gas is confined partly by changing the size of the 
confining potential in one degree of freedom. According 
to the theoretical developments of Sec. II it is assumed 
that the resulting physical effects of such a confine- 
ment can be described by periodic boundary conditions. 
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FIG. 1. Dependence of the critical temperature T c of an in- 
teracting Bose gas on the scattering length a. (d — (V/N) 1 ^ 3 
denotes the mean distance between atoms and T is the crit- 
ical temperature of an ideal Bose gas.) 

Let us first of all concentrate on the critical proper- 
ties of a homogeneous interacting Bose gas in three spa- 
tial dimensions. Though the corresponding characteristic 
critical exponents are already well known jl9| , non uni- 
versal thermodynamic properties of a three-dimensional 
homogeneous Bose gas are still subject to controversial 
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discussions |l^,|ri|,|l3| . One example of such a thermody- 
namic property which is of topical interest is the critical 
temperature and its dependence on the scattering length 
a characterizing the strength of the interparticle interac- 
tion. In three spatial dimensions this scattering length 
is related to the characteristic interaction parameter G 
of Sec. II by G = AnaK. The ratio between the crit- 
ical temperature T c of an interacting Bose gas and the 
corresponding value of an ideal Bose gase T can be de- 
termined easily from Eq. (p2[) . with the help of the re- 
lation T c /T = [(iV/V)A^]/2.612 (Thereby the relation 
(N/V)Xf h = C(3/2) w 2.612 has been used which applies 
to an ideal homogeneous Bose gas at the transition tem- 
perature.) The resulting predictions based on the RG 
approach of Sec. II are depicted in Fig.l. One realizes 
that in this particular case of a repulsive interparticle in- 
teraction, i.e. a > 0, the critical temperature increases 
with increasing scattering length. The quantitative de- 
pendence shown in Fig. 1 is approximated to a high degree 
of accuracy by the polynomial 

T c = T o [1.000 + 3.423(a/d) - 29.986(a/d) 2 + 

145.183(o/d) 3 ]. (25) 

From this approximation we conclude that in the 
limit of a vanishing interparticle interaction, i.e. 
a -> 0, one obtains [(T c - T )/T ] = 3.423(a/d) + 
0[(a/d) 2 ]. Eq.(p5|) is in satisfactory agreement with 
recent experimental results |l2|| and with the per- 
turbative and RG results of Refs. |0Jl| and @. 
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FIG. 2. Dependence of the scaled critical density 
[(N/V)\th\ °f an interacting Bose gas on the scattering length 
a. (Xth denotes the thermal wave length.) 

The critical density and its monotonically decreasing 



dependence on the scattering length is shown in Fig. 2. 
Note that in the limit a — > one obtains the result 
(N/V)Xf h — > 2.612 as expected for an ideal Bose gas. 

The critical quantities depicted in Figs.l and 2 can 
be obtained from the RG approach described in Sec. 
II in a straightforward way from the knowledge of the 
critical trajectory. This critical trajectory is the sta- 
ble manifold of the unstable fixed point for D = 3. 
This unstable fixed point is given by M c = 1/2 and 
g c = tt 2 /2 with g(l) = G(l)/(j3(l)/j3 A ). The associ- 
ated eigenvalues of the linearized RG-equations for M(l) 
and g(l) are A x = (3 - V33)/6 and A 2 = (3 + V33)/6. 
The latter eigenvalue which is associated with the unsta- 
ble manifold determines the various critical exponents 
of the second order phase transition. Thus the scaling 
of the correlations length of the interacting Bose gas, 
for example, is determined by the characteristic expo- 
nent v — I/A2 = 0.686139 which compares well with the 
known result of v = 0.67 [ [l9f which is based on a per- 
turbative approach to the RG above the transition point. 
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FIG. 3. Dependence of the second order coherence factor 
gr' 2 '(0) of an interacting Bose gas on the scattering length a. 

One of the quantities which is accessible experimen- 
tally and which has received attention recently j^] is the 
spatially averaged second order coherence factor of an 
interacting Bose gas. This quantity can be determined 
directly from the RG equations (15-17) of Sec. II by 
using relation (pif). Its critical value at the transition 
point and its dependence on the scattering length are 
depicted in Fig. 3. For vanishingly small values of the 
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interatomic interaction this second order coherence fac- 
tor approaches the value of two. This particular value is 
also known to characterize photon bunching of a chaotic 
electromagnetic field. With increasing scattering length 
a the critical value of the second order coherence factor 
decreases thus indicating that with increasing repulsive 
interactions the bosons tend to avoid each other. 

Let us now investigate the influence of confining this 
homogeneous interacting Bose gas with respect to one de- 
gree of freedom, say the z-direction. Furthermore, let us 
assume that the influence of this confinement can be de- 
scribed quantitatively by periodic boundary conditions. 
In such a case the resulting density of states d(l) which 
enters the RG equations (15-17) is given by Eq.(pO|). In 
Fig. 4 isotherms of the pressure of the interacting Bose 
gas are depicted for various values of the characteristic 
confinement length L z in the z-direction. Thereby above 
the phase transition, i.e. for negative values of the scaled 
chemical potential M(l), renormalization group equa- 
tions have been used which apply to a vanishing order 
parameter M. 
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FIG. 4. Isotherms of an interacting Bose gas for various 
values of the characteristic length of confinement L z . The 
scaled pressure [pA? h /(fcT)] is depicted as a function of the 
scaled volume [(V/7V)A t ~ 3 ]. 

The temperature of these isotherms is chosen so that 
Xth = lOV^ra ~ 25a. For Rubidium atoms, for exam- 
ple, with a scattering length of magnitude a = 5.3nm 
this condition corresponds to a temperature of T = 
1.98 x 10 _6 K. If the length of confinement L z is large 
in comparison with both the thermal wave length X t h 



and the scattering length a (compare with full curve), 
the characteristic signatures of a well pronounced second 
order phase transition are realized at a critical volume 
of mag nitude (y/N)\-£ = 0.456 in agreement with the 
result of Fig. 2. As soon as the length of confinement L z 
becomes comparable to the thermal wave length \ t h the 
pressure dependence is modified significantly. The most 
prominent feature of the depicted pressure dependence is 
the disappearance of the characteristic signature of the 
second order phase transition and the smoothing of this 
pressure dependence. 

In view of the recent interest in the behavior of 
the second order coherence factors of interacting Bose 
gases its isothermal dependence on the density of 
the interacting Bose gas is depicted in Fig. 5. 
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FIG. 5. Isotherms of the second order coherence factor 
<T ^(0) of an interacting Bose gas. 

The parameters chosen are the same as the ones used in 
Fig. 4 . Again for very large values of the confinement 
length L z one notices a well pronounced second order 
phase transition at the same critical volume as in Fig. 4. 
Far below the transition point, i.e. for very small densi- 
ties, the second order coherence factor approaches a value 
of one which is the expected value for a Bose condensate 
at zero temperature. Lowering the length of confinement 
L z in the z-direction the characteristic signature of the 
second order phase transition disappears. 
sectionSummary 

In summary, a theoretical description of thermody- 
namic properties of an interacting confined Bose gas 
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has been developed which is based on a RG approach. 
Thereby physical effects which originate from the pres- 
ence of confinement have been taken into account by pe- 
riodic boundary conditions. Thus this approach yields a 
proper description of thermodynamic properties also in 
cases in which the characteristic length of confinement 
becomes comparable to the thermal wave length of a 
Bose gas and in which local density approximations break 
down. A more realistic and more complete treatment 
of boundary conditions will be postponed to subsequent 
work. It has been demonstrated that this approach gives 
a reliable description of the critical behavior of an inter- 
acting Bose gas which is consistent with recent theoreti- 
cal and experimental work. Furthermore, concentrating 
on the bunching properties of an interacting Bose gas the 
influence of the confinement length on the characteristic 
aspects of the phase transition has been worked out. 

This work is supported by the Deutsche Forschungsge- 
mcinschaft within the Forschergruppe 'Quantengase'. 



5(4, 4*) = 5(4<, 4* ) + ^5<p T M5<p (27) 

with 

Sip T = (...<V m „, (Sip mn )*,Sip 

— m—n )*,...). (28) 

Thereby one chooses the symmetry breaking field Tp iden- 
tical to the most probable uniform field configuration of 
Eq.(|l2|) so that terms linear in vanish in Eq. (p6|) . 
The small-momentum field components tp m n which are 
not integrated out and which are kept constant consti- 
tute the field 4<. The symmetric matrix M is given by 



M 



with 
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A — iuj 
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A — iLO n 






(29) 



APPENDIX A 

In this appendix the derivation of the RG transfor- 
mation of Sec. II is outlined in which an infinitesimal 
momentum shell of the partition function Z(z, (3) is inte- 
grated out. These equations are valid in the dynamical 
regime below the transition temperature which is char- 
acterized by a non-zero order parameter and a positive 
renormalized chemical potential M(l). 

Starting point is the path integral representation of 
Eq. (U) with the (scaled) action functional of Eq. (g) . In- 
troducing a symmetry breaking, real-valued uniform field 
configuration Tp this action functional can be written in 
the form 

5(4, 4*) = {-Mip 2 + \gTp±} + [^00 + ip* 00 ][-M7p + gep 3 } + 



m,n 



M + 2gTp 2 ]}\<p 



iii. 



mi,ni...m 3 ,n3 



< 

e> - M + 2g[Tp 2 +Tp^ 



mn 



< < 



B = g[<f 2 



< 



mn 



< 

E 

mn 



(30) 



The quantity e> = 1/2 denotes the scaled energy of the 
eliminated momentum shell. According to the RPA [fl4|| , 
in Eq.(p7|) only products of large-momentum field com- 
ponents Sip mn are kept for which the sum of the momenta 
and of the Matsubara frequencies vanishes. Performing 
the Gaussian integrations over the fields 5tp mn one ob- 
tains 



Z(z,{3) = [H J 



N n 



mn i — ,s 
e ™ 



(31) 



<i5l mi _„ 1 ^I m2 _„ 2 -^ m 3„ 3 ](5(n 1 +n 2 + n 3 )5(m 1 + m 2 + m 3 ) 
1 
2 



with the new effective action 

5 ne 44<,4*) = 5(4<,4*)- 

(Cm„/V/3) 2 + U 2 n 



It, 



III 71 — — OQ 



ln- 



A 2 



E 99 
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^m 4 n 4 ^(ni + n 2 - 



X 



-mi -ni ¥-m 2 -ni <Pm 3 n 3 

- n 3 + n 4 )(5(mi + m 2 + m 3 + 1114) 



5(4<,4*<)- 



1 



^ln 6 



u>l-\B\* 

(/3//3a)(e>-M) 



with g = (VA D )~ 1 (P/(3\)~ 1 G. In order to eliminate 
the infinitesimal momentum shell around the maximum 
momentum (SA), for which e~ l < |k m | < 1 (0 < I <C 1), 
one expands the (scaled) action functional of Eq.@ up 
to second order in the large-momentum field components 
StPmn, i-e. 



(26) 

Thereby the relation 
sinh 2 (x)/a; 2 



2^ 4sinh 2 [A(4<)/2]- 



(32) 



n 



n— — oo,n^0 



+ ) 2 ] 

nir 



has been used to perform the summation over all Mat- 
subara frequencies. The quantity A(4<) is given by 
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A($<) = (/?//3a) Ve>[e> + 2Af] - A($<) (33) 



with 



mn / J ~r 

m n— — oo 

< oo 

■Py^ E [VmnV-m-n + (¥>mn)*(</>-m-n)*] 

111 71 — — OO 

< OO 

(2£> + C/v5)E X! Iv« 



- mn | ? 



m n — — oo 



C = (-4e> - 2M)M/Tp, 
D = -3M 2 /lp 2 . 



(34) 



The symbol < in the summations of Eq.([54h indicates 
that only momenta have to be taken into account which 
have not yet been integrated out. The prime in the sum 
of Eq.(|32|) indicates summation over all the eliminated 
momentum components. In the continuum limit this lat- 
ter summation reduces to 



E = ( VA °) 



(2^ 



-I (0<i<l). 



(35) 



For the derivation of the RG transformation we first of 
all split off the zero-temperature contribution according 
to 



8S(M, (3). 



(/3//3 A )«n(M,/3-oo) + 
(36) 



The remaining temperature dependent contribution 
5S(M,(3) is then expanded up to second order in the 
small-momentum field components which constitute the 
field <&<, i.e. 



5S(M, (3) 



2 4sinli (A/2) 

1 Q3//? A ) 2 A($<) ' d 

2 2A ^ ' 

m 

l r (/3//3 A ) 2 A($ 



2A 



d\ 4sinh 2 (A/2) 
Id 2 ~ x 



A ^ = m E(' 3 // ? A) 2 ^[ coth ( A / 2 ) - w 



2A 1 



(39) 



Eq.(p7|) implies a change of the chemical potential, i.e. 



(40) 



with 



AM = - V(/3//3 A ) 2 7?r[coth(A/2) - l]D + 2^-M - 
L — ' 2A if 



1 



~[coth(A/2) - 1]}C 2 . (41) 
A 



The quantities C and D are defined by Eq. (p4j) . The 
corresponding change of the scaled two-body coupling 
strength G is determined by the requirement that 



V new = yJ(VAV)([3/(3 A )M new /G, n 



(42) 



in accordance with Eq. ( |l2| ) . This implies that 

G new = G + AG (43) 



with 



(44) 



Eqs.(p8| - and (|3| - 0) characterize the Kadanoff 
transformation [|J, i.e. the elimination of the infinites- 
imal momentum shell for which e~ l < |k m | < 1 with 
( < 1. The RG transformation is completed by perform- 
ing the scaling transformations which restore the original 
momentum cut-off (hA), namely 

km' = k m e , 

V = Ve~ Dl , 
= Pe- 21 , 



(45) 



2dA 2 4sinh 2 (A/2) 



with A = A($< = 0) and with A($<) as defined by 
Eq. (pi[) . Note that the terms involving e A result from 
the separation of the zero-temperature contribution ac- 
cording to Eq. (|T3[) . This expanded effective action in- 
volves terms linear in the field <&<. These linear terms 
can be absorbed in a redefinition of the most probable 
configuration, i.e. 



With M(l) = (Af + AAf)e 2 ', G(l) = (G + AG)e< 2 - D ", 
{3? f) (I) — f3e~ 21 these scaling relations together with the 
Kadanoff transformations of Eqs.(^8|- |d],[l3]- ^J) imply 
the RG equations of Eqs.([l5| - pi). 



Vnew =<P + &f 



with 



APPENDIX B 

For the sake of completeness in this appendix the RG 
equations are summarized which apply to the case of a 
vanishing order parameter above the transition point. 
(38) In this dynamical regime the scaled chemical potential 

M(l) might become negative. The RG equations which 
apply to this regime of negative values of M(l) can be 
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obtained in a straightforward way with the help of second [15] R. P. Feynman, Statistical Mechanics, (Benjamin, Read- 
order perturbation theory by assuming a vanishing order ing, 1972). 

parameter §,0. Thus, for M(l) < one obtains the RG [16] F. J. Wegner and A. Houghton, Phys.Rev. 8, 401 (1973). 
flow equations [17] F. W. Wiegel and J. B. Jalickee, Physica 57, 317 (1972). 

[18] A. L. Fetter and J. D. Walecka, Quantum Theory of Many 
du)(l) — Dh ( -(/3//3a)(c -M{1))\ Particle Systems, Me Graw-Hill, New York, 1971; chs. 6, 

~dT = WM > 10 - 

dM(l) 1 ^ Zinn- Justin, Quantum Field Theory and Critical Phe- 

= 2M(l) - d(l)G(l){ 1}, nomena, Oxford University Press, New York, 1989. 



dl ww ^ y L tanh[(/3/ / 3 A )(e > - M(l))/2] 

dG{l) - (2-D)G(l)-d(l)G 2 {l) x 



dl 

1 1 

i tanh[(/3//3 A )(e > - M(/))/2] 2(/?//3 A ) 2 (e> - M{1)) + 

(/?/%) 5 }. (46) 

Starting from a positive chemical potential, i.e. M(l — 
0) = M > 0, one has to switch from Eqs. (16-17) to Eq. 
( [46| ) as soon as M{1) becomes negative in the process of 
renormalization. If one starts from a negative chemical 
potential, i.e. M (I = 0) = M < 0, one has to solve Eq. 
( (4(f ) with the intial conditions 

u(l = 0) =0, 
M(l = 0) = M < 0, 
G(l = 0) = G. (47) 

The corresponding grand thermodynamic potential is 
given by Q(M, 0) = uj{1 -> oo). 
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